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ABSTRACT 

We study the use of red sequence selected galaxy spectroscopy for unbiased estimation of galaxy clus- 
ter masses. We use the publicly available galaxy ca talog produced usin g the semi-analytic model of 
iDe Lucia & Blaizoj (l2QQ7h on the Millenium Simulation dSpringel et alJ 120051) . We make mock observations 
to mimic the selection of the galaxy sample, the interloper rejection and the dispersion measurements for large 
numbers of simulated clusters spanning a wide range in mass and redshift. We explore the impacts on selection 
using galaxy color, projected separation from the cluster center, and galaxy luminosity. We probe for biases 
and characterize sources of scatter in the relationship between cluster virial mass and velocity dispersion. We 
identify and characterize the following sources of bias and scatter: intrinsic properties of halos in the form of 
halo triaxiality, dynamical friction of red luminous galaxies and interlopers. We show that due to halo triaxial- 
ity the intrinsic scatter of estimated line-of-sight dynamical mass is about three times larger (30 — 40%) than 
the one estimated using the 3D velocity dispersion (~ 12%) and a small bias (< 1%) is induced. Furthermore 
we find evidence of increasing scatter as a function of redshift and provide a fitting formula to account for 
it. We characterize the amount of bias and scatter introduced by dynamical friction when using subsamples of 
red-luminous galaxies to estimate the velocity dispersion. We study the presence of interlopers in spectroscopic 
samples and their effect on the estimated cluster dynamical mass. Our results show that while cluster velocity 
dispersions extracted from a few dozen red sequence selected galaxies do not provide precise masses on a sin- 
gle cluster basis, an ensemble of cluster velocity dispersions can be combined to produce a precise calibration 
of a cluster survey mass-observable relation. Currently, disagreements in the literature on simulated subhalo 
velocity dispersion- mass relations place a systematic floor on velocity dispersion mass calibration at the 15% 
level in mass. We show that the selection related uncertainties are small by comparison, providing hope that 
with further improvements to numerical studies this systematic floor can be substantially reduced. 
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1. INTRODUCTION 

Clusters of galaxies are the most massive collapsed objects 
in the Universe and sit at the top of the hierarchy of non- 
linear structures. They were first identified as over-den se re- 
sions in the projected number counts of galaxies (e.g. Abelll 
[1958, Zwicky et al. 1968). However nowaday s clusters can be 
identified over the whole electro-magnetic range, including as 
X-ray sources (e.g. .Bohringer et al. 2()00[ iPacaud et al.ll2007L 
IVikhlininet"ani2009L ISuhada et al. 2012), as o ptically over- 
densi ti es of red galaxies ('Gladders & Yee 2005llKoester et al.l 
I2007L iHao et a l. 2010, Szaboetal. 20lll) and as distor- 
sions of the cosmic ni icrowave backgrou nd as predicted by 
ISunvaev & Zel'dovichl (197 2} (e.g. Vanderlindee t al.ll2010L 
[Marriage et al. 201 1, Planck Collaboration et al. 20TT|). 

Given the continuous improvement in both spatial and spec- 
tral resolution power of modern X-ray, optical and infrared 
telescopes, more and more details on the inner properties of 
galaxy clusters have been unveiled in the last decade. These 
objects, that in a first approximation were thought to be viri- 
alized and spherically symmetric, have very complex dy- 
namical features - such as stro ng asymmetries and clumpi- 
ness (e.g. Geller & Beerslll982l Dressier & Shectmanlll988i 
iMohr et al . 1995) - witnessing for violent processes being 
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acting or having just played a role. They exhibit luminos- 
ity and temperature functions which are not trivially related 
to their mass function, as one would expect for virialized 
gravitation-driven objects. Moreover, the radial structure of 
baryons' properties is far from being completely understood: 
a number of observational facts pose a real challenge to our 
ability in modeling the physics of the intracluster medium and 
the closely related physics of the galaxy population. Indeed 
a number of different physical processes are acting together 
during the formation and evolution of galaxy clusters. Gas 
cooling, star formation, chemical enrichment, feedback from 
supemovae explosions and from active galactic nuclei, etc are 
physical processes at the base of galaxy formation, which are 
difficult to disentangle (e.g. see lBensonll201 for a recent re- 
view on galaxy formation models). 

Line-of-sight galaxy velocities in principle provide a mea- 
sure of the depth of the gravitational potential well and there- 
fore can be used to estimate cluster masses. Furthermore, 
galaxy dynamics are expected to be less affected by the com- 
plex baryonic physics affecting the intra cluster medium. 
Thus, one would naively expect a mass function defined on 
the basis of velocity dispersion to be a good proxy of the un- 
derlying cluster mass. However a number of possible sys- 
tematics can affect dynamic al mass estimation a nd must be 
carefully take into account. iBiviano et al.l (l2006h for exam- 
ple studied a sample of 62 clusters at redshift z = from 
a ACDM cosmological hydrodynamical simulation. They 
estimated virial masses from both dark matter (DM) particles 
and simulated galaxies in two independent ways: a virial mass 
estimator corrected for the surface pressure term, and a mas s 
estimator based entirely on the velocity dispersion ay . They 
also modeled interlopers by selecting galaxies within cylin- 
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ders of different radius and length 192/i~^Mpc and applying 
interloper removal techniques. They found that the mass es- 
timator based entirely on velocity dispersions is less sensitive 
on the radially dependent incompleteness. Furthermore the 
effect of interlopers is smaller if early type galaxies, defined 
in the simulations using their mean redshift of formation, are 
selected. However, the velocity dispersion of early t ype galax- 
ies is biased low with respect to DM particles. E vrard et al.l 
(l2008h analysed a set of different simulations with different 
cosmologies, physics and resolutions and found that the 3D 
velocity dispersion of DM particles within the virial radius 
can be expressed as a tight function of the halo virial masfl 
regardless of the simulation details. They also found the scat- 
ter about the mean relation is nearly log-normal with a low 
standard d eviation crino- — 0.04. In a more recent work, 
I White et a l. (2010) used high resolution N-body simulations 
to study how the influence of large scale structure could affect 
different physical probes, including the velocity dispersion 
based upon sub-halo dynamics. They found that the highly 
anisotropic nature of infall material into clusters of galax- 
ies and their intrinsic triaxiality is responsible for the large 
variance of the ID velocity dispersion under different lines 
of sight. They also studied how different interloper removal 
techniques affect the velocity dispersion and the stability of 
velocity dispersion as a function of the number of sub-halos 
used to estimate it. They found that only when using small 
numbers of sub-halos (< 30) is the line of sight velocity dis- 
persion biased low and the scatter significantly increases with 
respect to the DM velocity dispersion. Furthermore the ef- 
fect of interlopers is different for different interloper rejection 
techniques and can significantly increase the scatter and bias 
low velocity dispersion estimates. 

Currently IR, SZE and X-ray cluster surveys are de- 
livering significa nt numbers of clusters at redshifts 

z > 1 (e.g. [Stanford et a 

20Q9L iFa ssbende r et al.1 12011 
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[Williamson etaH |2QT1 
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samples is challenging using weak lensing, making velocity 
dispersion mass estimates particularly valuable. At these red- 
shifts it is also prohibitively expensive to obtain spectroscopy 
of large samples of cluster galaxies, and therefore dispersion 
measurements must rely on small samples of 20 to 30 cluster 
members. This makes it critically important to understand 
how one can best use the dynamical information of a few 
dozen of the most luminous cluster galaxies to constrain the 
cluster mass. It is clear that with such a small sample one 
cannot obtain precise mass estimates of individual clusters. 
However, for mass calibration of a cluster SZE survey, for 
example, an unbiased mass estimator with a large statistical 
uncertainty is still valuable. 

In this work we focus on the characterisation of dynamical 
mass of clusters with particular emphasis on high-z clusters 
with a small number of measured galaxy velocities. The plan 
of the paper is as follows. In Sec. [2] we briefly introduce the 
simulation describe the adopted semi-analytic model, and in 
Sec. [3] we present the results of our analysis. Finally, in Sec. 
[51 we summarise our findings and give our conclusions. 

2. INPUT SIMULATION 



^ Throughout the text, we will refer to Mvir as the mass contained within 
a radius i?vir encompassing a mean density equal to 200 pc, where pc is the 



TABLE 1 

The redshift-number distribution of the 22,484 clusters with 

mvir > lo^'^ m0 analysed in this work at different redshift. 

Column 1: redshift z; Column 2: number of clusters A^cius- 
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critical cosmic density. 



This analysis is based on the publicly available galaxy 
cat alogue produced using th e semi-analytic model (SAM) 
by De Lucia & Blaizo^ (I2QQ7I) on the Millennium Simulation 
( Springel et al. 2005). The Millennium Simulation adopts the 
following values for the parameters of a flat A cold dark mat- 
ter model: VLdm = 0.205 and Q^ = 0.045 for the den- 
sities in cold dark matter and baryons at redshift z = 0, 
(78 = 0.9 for the rms linear mass fluctuation in a sphere of ra- 
dius Sh~^Mpc , h = 0.73 for the present dimensionless value 
of the Hubble constant and n = 1 for the spectral index of the 
primordial fluctuation. The simulation follows the evolution 
of 2160^ dark matter particles from z = 127 to the present day 
within a cubic box of 500/i~^Mpc on a side. The individual 
dark matter particle mass is 8.6 x 10^h~^ Mq. The simula- 
tion was carried out with the massively parallel GADGET-2 
code (Springel 2005). Gravitational forces were computed 
with the TreePM method, where long-range forces are cal- 
culated with a classical particle-mesh method while short- 
range forces are dete rmined with a hierarchical tree approach 
dBames & H ut!' 1986"). The gravitational force has a Plummer- 
equivalent comoving softening of 5/i~^kpc, which can be 
taken as the spatial resolution of the simulation. Full data 
are stored 64 times spaced approximately equally in the loga- 
rithm of the expansion factor. Dark matter halos a nd subhalos 
were identified with th e friends-of- friends (F OF; 'Davis et aH 
IT985h ) and SUBFIND dSpringel et al.|[2QQlah algorithms, re- 
spectively. Based on the halos and subhalos within all the 
simulation outputs, detailed merger history trees were con- 
structed, which form the basic input required by subsequently 
applied semi-analytic models of galaxy formation. 

We recall that the SAM we emplo y builds upon the n iethod - 
ology originally introduced by Kauffmann et al. ('1999^, 
Springel et al. (2001b) and De Lucia et al. (2004b). We refer 
to the original papers for details. 

The SAM adopted in this study includes explicitly DM sub- 
structures. This means that the halos within which galaxies 
form are still followed even when accreted onto larger sys- 
tems. As explained in Springel et al. (2001) and De Lucia 
et al. (2004), the adoption of this particular scheme leads 
to the definition of different galaxy types. Each FOF group 
hosts a Central galaxy. This galaxy is located at the posi- 
tion of the most bound particle of the main halo, and it is 
the only galaxy fed by radiative cooling from the surrounding 
hot halo medium. Besides central galaxies, all galaxies at- 
tached to DM substructures are considered as satellite galax- 
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Fig. 1 . — The evolution of the normahzation A (left panel) and slope B (right panel) parameters used to fit the relation between the 3D velocity dispersion of 
all the galaxies within i^vir^nd the virial mass of each cluster (eq. |2). Red horizontal dashed lines represent the mean value. The slope is moderately shallower 
than the self-similar expectation. 



ies. These galaxies were previously central galaxies of a halo 
that merged to form the larger system in which they currently 
reside. The positions and velocities of these galaxies are fol- 
lowed by tracing the surviving core of the parent halo. The 
hot reservoir originally associated with the galaxy is assumed 
to be kinematically stripped at the time of accretion and is 
added to the hot component of the new main halo. Tidal trun- 
cation and stripping rapidly reduce the mass of DM substruc- 
tures (but not the associated stellar mass) below the resolu- 
tion limit of the simulation ( De Lucia et al.|[2004al: Gao et al. 
2004). When this happens, we estimate a residual surviving 
time for the satellite galaxies using the classical dynamical 
friction formula, and we follow the positions and velocities 
of the galaxies by tracing the most bound particles of the de- 
stroyed substructures. 

3. PROPERTIES OF THE FULL GALAXY POPULATION 

3.1. Intrinsic galaxy velocity dispersion 

lEvrard et"an (1200 81) showed that massive dark matter halos 
adhere to a virial scaling relation when one expresses the ve- 
locity dispersion of the DM particles as a function of the virial 
mass of the halo in the form: 



0'DM{Myir^z) = CFDM.lb 



h{z)Myir 
1015 M0 



(1) 



where aDM,i5 = 1082.9 ± 4.0 km s ^ is the typical 3D 
velocity dispersion of the DM particles within i^vir for a 
10^^ h-^ Mo cluster at z = and a = 0.3361 ±0.0026. Sim- 
ilarly, we first compute for each cluster the 3D velocity disper- 
sion asD (divided by >/3) of all the galaxies within i^virand 
then fit the relation between asD and Mvir in the form of 
log{crsD) oc /o^(/i7o(^)Mvir/lO^5M0) individually at any 
of the redshift listed in Table[T] As a result we can express the 
dynamical mass Mdyn as: 



Mdyn= [^Y h7o{z)-Ho''MQ, 



(2) 



where the resulting best fitting values of A and B with their 
associated error-bars are shown in Fig. [Has a function of 



redshift. Dashed horizontal red lines show the average jval- 
ues which are respectively A = 938 ± 3kms~i and B = 
2.91 ±0.03. 

After accounting for the differences in the Hubble pa- 
rameter, our measured normalization of the g alaxy veloc- 
ity dis persion- mass relation is within ^3% of lEvrardet al.l 
(|2008|) . This reflects the differences between the subhalo and 
DM particle dynamics. As has been previously pointed out 
(e.g. taao et al. 2 004, Goto 2005.Faltenb acher & Diemanl 
2006."Evra rd et al. 11200 SllWhite et al.ll2010h . the velocity bias 
between galaxies and DM is expected to be small by < 5%. 
But to be absolutely clear, we adopt our measured galaxy ve- 
locity dispersion- mass calibration in the analyses that fol- 
low. To better visualize the relative importance of the cos- 
mological redshift dependence we show in Fig. [2] the redshift 
evolution of the normalisation parameter A (solid black line 
)when the fit is made on the relation log{(TsD) o^ log{h7o{z = 
0)Mvir /lO^^M©). The expected self- similar evolution given 
hy A{z) = Ax E{z)^ is highlighted (dashed red line), where 
the term A is equal to mean value A = 938 km s"^ and E{z) 
describes the universal expansion history H{z) = HqE{z). 
In other words. Fig. [2] shows the typical galaxy velocity dis- 
persion in kms~i for a cluster with Mvir = 10^^ Hj^Mq as 
a function of redshift and demonstrates the nearly self-similar 
evolution (within ~1%) over the redshift range tested in this 
work. 

For the full sample of clusters analysed (see Table [TJ, we 
then compute the dynamical masses by applying Eq.[2lto (1) 
the 3D galaxy velocity dispersion (divided by >/3) and (2) 
to each orthogonal projected ID velocity dispersion. Fig. [3] 
shows the comparison between the virial masses Mvir and 
the resulting dynamical masses M^d (left panel) and Mid 
(right panel) for the full sample of clusters. The best fit of 
the relation (dashed black and white lines) is virtually indis- 
tinguishable from the one-to-one relation (dotted-dashed pur- 
ple line) in the case of the 3D velocity dispersion. On the 
other hand, in the case of the ID velocity dispersion there is 
a small but detectable difference between the one-to-one re- 
lation and the best fit. The best fit of the dynamical mass 
for the ID velocity dispersion is about < 1% lower than the 
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Fig. 2. — The evolution of the normalisation parameter A{^z) of Eq.[2lwhen 
no self-similar evolution is taken into account (solid black line). Red dashed 

line is showing the best fitting parameter A x {H{z)/H{))^ 



one-to-one relation. We will show in Section [3^ that this dif- 
ference can be explained in terms of triaxial properties of ha- 
los. Typical logarithmic scatter of (Ju-iDlM^i^ — 0.145 and 
cTMiD/M^ir. — 0.334 are highlighted with dotted black and 
white lines in Ioqiq scale. We find that, similar to results by 
IWhite et a l. (2010), using the ID velocity dispersion rather 
than the 3D velocity dispersion increases the intrinsic log scat- 
ter around the mean relation by a factor of '^ 2.3. 

We further investigate the intrinsic scatter in the relation 
between the true virial masses and the dynamical mass esti- 
mates in Fig. m Taking a to be the standard deviation of 
the logarithm of the ratio between the dynamical mass esti- 
mate and the virial masses, we show that in the case of the 3D 
velocity dispersion (dashed red line) and the ID velocity dis- 
persion (dotted black line) the scatter increases with redshift. 
The solid black line shows a linear fit to the evolution of the 
intrinsic Mdyn-iD scatter and can be expressed as: 



^ln{MiD/M^,r-) 



0.3 + 0.075 X z. 



(3) 



Velocity dispersions are '^25% less accurate for estimating 
single cluster masses at z = 1 than at low redshift. 

The logarithmic scatter of the ID velocity dispersion mass 
estimator a Mm around the true mass arises from two sources 
of scatter: (1) the logarithmic scatter between the 3D velocity 
dispersion mass estimator and the true mass - a[MsD/Myir] 
(red dashed line in Fig |4]) and (2) the logarithmic scatter be- 
tween the ID and 3D velocity dispersions crlaiD/osD] (solid 
green line). The expected ID dispersion mass scatter is then 
the quadrature addition of these two sources: 



(^M,n ^ (^^[MSD/Myir] + {5 X aiaiD/cTSDJV 



(4) 



where B is the best fitting slope parameter from Eq.|2] The ex- 
pected a Mid estimate from Eqn.|4] appears as a dotted-dashed 
purple line in Fig [H note that this estimate is in excellent 
agreement with the directly measured scatter ( dotted black 
line). Therefore, we show- as pointed out by 'Whit e et aH 
f2010)- that the dominant contributor to the scatter is the in- 
trinsic triaxial structure of halos. Furthermore its evolution 
with redshift is also the dominant source of the increasing 
scatter of the ID dynamical mass estimates with redshift. By 
comparison, the scatter between the 3D velocity dispersion 
mass estimator and the true mass a[ln{MsD/Myir)], which is 



reflecting departures from dynamical equilibrium due to on- 
going merging in the cluster population, is relatively minor. 
Ultimately it is the lack of observational access to the full 3D 
dynamics and distribution of the galaxies that limits us from 
precise single cluster dynamical mass estimates. 

3.2. Triaxiality 

The presence of pronounced departures from sphericit y in 
dark matter halos (Thomas & Couchman 1992, Warren e t al.l 
II992L Ijing & Suto 2002), if not approximately balanced be- 
tween prolate and oblate systems, could in principle not only 
increase the scatter in dynamical mass estimates, but also lead 
to a bias. If, for example, clusters were mainly prolate sys- 
tems, with one major axis associated to a higher velocity dis- 
persion and two minor axes with a lower velocity dispersion, 
there should be two lines of sight over three associated with 
a lower velocity dispersion. This could potentially lead to a 
bias in the ID velocity dispersion with respect to the 3D ve- 
locity dispersion. To quantify this possible bias, we compute 
the moment of inertia for each cluster in the sample, and we 
then calculate the velocity dispersions along each of the three 
major axes. As has been pointed o ut before (lTormenlll997i 
iKasun & Evr^l2QQ5L IWhite et al.1 .2010) the inertia and ve- 
locity tensor are quite well aligned, with typical misalignment 
angle of less than 30°. In Fig. [51 at each redshift we show the 
lowest velocity dispersion cfq with black crosses, the highest 
(72 with green stars and the intermediate one ai with red dia- 
monds normalized to the 3D velocity dispersion cj^d (divided 
by v^S). Dashed blue lines are the 16, 50 and 84 percentile of 
the full distribution and DEV is the associated standard devia- 
tion which, as expected from Fig.|4]is increasing with redshift. 
A perfectly spherical cluster in this plot will therefore appear 
with the three points lying all at the value 1 , whereas prolate 
and oblate systems will have the intermediate velocity disper- 
sions (Ji closer to the lower one cjq and to the higher one (72, 
respectively. The black solid line is the best fit of the distribu- 
tion of the intermediate ai velocity dispersions and it is very 
close to unity, showing that dynamically, clusters do not have 
a very strong preference among prolate and oblate systems. 
Furthermore this result is true for the range of redshifts and 
masses we examine here. 

This can be better seen in Fig. [51 where we show that we 
measure only a mild excess (at < 5% level) of prolate sys- 
tems. In addition, for each cluster in the sample, we compute 
a "prolateness" quantity Prol as: 



Prol 



((72 -CFi)- ((71 -(7o) 



cr^D 



(5) 



A prolate system will have a positive Prol value whereas an 
oblate one will have a negative Prol. Fig. [6] shows a map 
representing the distribution of the Prol variable as a function 
of the cluster mass (left panel) and redshift (right panel). To 
compute the former we stack clusters from all the redshifts, 
and for the latter we stack clusters from all masses. As it is 
shown in Fig. [6l there are no clear dependencies of the Prol 
variable on the cluster mass or redshift. The slight excess of 
prolate over oblate systems at all masses and redshifts would 
translate into ID dynamical masses slightly biased towards 
smaller masses. Indeed, this is seen as a ~ 1% effect in Fig.[3l 

4. PROPERTIES OF SPECTROSCOPIC SUBSAMPLES 

Results in the previous section relied on the full galaxy sam- 
ple within each cluster virial region. We now study the pos- 
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Fig. 3. — The relation between Mvir and the dynamical mass for all the clusters in the sample. For each cluster the dynamical mass is inferred by applying Eq. 
[21 to the 3D velocity dispersion divided by a/3 (left panel) and for each of the three projected ID velocity dispersions (right panel) of all the galaxies extracted 
from the De Lucia & Blaizot (2007) database within i?vir from the centre of the cluster. The dashed (dotted) white-black line is the best fit of the relation (plus 
and minus one a) and is virtually indistinguishable from the one-to-one relation (dotted-dashed purple line). 
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Fig. 4. — The redshift evolution of the logarithmic la scatter for the fol- 
lowing quantities: (1) the 3D galaxy velocity dispersion mass estimate scat- 
ter (dashed red), (2) the ID galaxy velocity dispersion mass estimate scatter 
(dotted black), (3) a fit to #2 (soHd black; Eqn.O, (4) the scatter of the ID 
velocity dispersion about the 3D dispersion (solid green), (5) the same quan- 
tity turned into mass scatter using Eqn. |2] (dashed-dotted blue) and (6) the 
expected ID dispersion mass scatter (#2) obtained by quadrature addition of 
#1 and #5, as explained in Sec. 13. IK dotted-dashed purple; Eqn.|4j. 



sible systematics affecting the cluster velocity dispersion and 
associated dynamical mass estimates when more realistic se- 
lection for the member galaxies are taken into account. 

We model the selection carried out in real world circum- 
stances by following the procedure we have developed for the 
South Pole Telescope (S PT) dynamical mass calibration pro- 
gram (iBazinet al.ll2Q12D . Namely, we preferentially choose 
the most luminous red sequence galaxies that lie projected 
near the center of the cluster for our spectroscopic sample. 
To do this we select galaxies according to their colors, which 



are a direct prediction of the adopted semi-analytic model. 
In particular, we compute the following color-magnitude di- 
agrams for different redshift ranges: ^ — r as a function of 
r for redshift z< 0.35, r — i as a function of i for redshifts 
0.35 <z< 0.75 and i — z sls sl function of z for redshifts 
larger than 0.75 (e.g. ISong et an i2011). We report in Fig. |7] 
the color-magnitude diagram at different redshifts for all the 
galaxies within the virial radius of each cluster. The model 
given by 'Song et al.l (12011 ). which has proven to be a good fit 
to the observational data, is highlighted with a dashed black- 
red line. As it is shown, the simulated cluster galaxy popula- 
tion has a red-sequence flatter than the observational results. 
Because the purpose of this work is not to study the evolution 
of the cluster galaxy population, but rather to see the effect of 
the selection of galaxies on the estimated dynamical mass, we 
adopt the following procedure: First we fit the red sequence at 
each analysed redshift. Then, we symmetrically increase the 
area on color-magnitude space in order to encompass 68% of 
the galaxies and iterate the fit. The resulting best fit and corre- 
sponding area are highlighted as green continuous and dashed 
lines in Fig. [T] Table [2l describes the width in color space used 
to select red sequence galaxies at each analysed redshift. 

This color selection helps to reduce the interlopers in our 
cluster spectroscopic sample. In addition to color selection, 
we explore the impact of imposing a maximum projected sep- 
aration R± from the cluster center, and we explore varying the 
spectroscopic sample size. In all cases we use the Ngai most 
massive (and therefore luminous) galaxies in our selected 
sample. Table [3] shows the range of Ngai and a = R^/vyir 
that we explore as well as the sample binning in redshift 
and mass. Note that for SZE selected clusters from SPT 
or equivalently X-ray selected samples of clusters, once one 
has the cluster photometric redshift one also has an estimate 
of the cluster virial mass and v irial radius from SZE signa- 
ture or X-ray luminosity (e.g. Reiprich & Bohringerl l2002[ 
Andersson et al. 201 1): therefore, we do in fact restrict our 
spectroscopic sample when building masks according to pro- 
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Fig. 5. — We show for eacli cluster tlie velocity dispersion along the tree major axes of the inertia momentum (black crosses for the smaller, red diamonds for 
the intermediate and green stars for the larger) normalized to the 3D velocity dispersion divided by \/3 as a function of the cluster mass in different redshift bins. 
The black solid line is the best fit of the intermediate axis velocity dispersion, and the dashed blue lines are the median and the 16 and 84 percentile of the full 
distribution. DEV is the associated standard deviation. 
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Fig. 6. — The distribution of the prolateness variable Prol (see Eqn.|5) as a function of the cluster mass for all the clusters at different redshift stacked together 
(left panel) and as a function of redshift (right panel). To guide the eye the dashed black and white line is highlighting the value of Prol = 0, while the dotted 
red- white lines are respectively the 16, 50 and 84 percentiles for the two different distributions. The cluster ensemble exhibits a slight preference for prolateness 
at all masses and redsliifts. 
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TABLE 2 

Color width of the simulated red sequence in magnitudes at 

each analysed redshift. column 1: redshift; column 2: 1 cr 

width of the red sequence in magnitudes. 

mag 



0.00 


0.05 
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0.06 


0.21 


0.08 


0.32 


0.10 


0.41 


0.06 


0.51 


0.08 


0.62 


0.08 


0.75 


0.08 


0.83 


0.09 


0.91 


0.09 


0.99 


0.07 


1.08 


0.06 
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0.05 



TABLE 3 

Parameter space explored for the mock observations. 

Column 1: maximum projected distance i?^ from cluster 

CENTER a = R^/Vvir]', COLUMN 2: Ng^l INITIAL NUMBER OF 
SELECTED MOST MASSIVE RED-SEQUENCE GALAXIES; COLUMN 3: 

REDSHIFT z; Column 4: cluster mass Mvir [10^^ M^; 
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jected distance from the cluster center relative to the cluster 
virial radius estimate. 

4.1. Dynamical friction and Velocity Bias 

In section 13.11 we showed the presence of a tight relation 
between the 3D dynamical mass and the virial mass Mvir 
for galaxy clusters. When dynamical masses are computed 
from the ID velocity dispersion instead of the 3D one, we 
significantly increase the scatter of this relation and intro- 
duce a negligible bias (< 1%) due to the triaxial properties 
of dark matter halos. We now study the effect of velocity 
segregation due to dynamical friction and its effect on the es- 
timated dynamical masses. To do this, for each cluster we 
select a number of red- sequence galaxies within the virial ra- 
dius i^virthat ranges from 10 to 100 galaxies as described in 
Table[3] We sort galaxies according to their luminosity (differ- 
ent bands were used at different redshift as described in Sec. 
|4]). This results in a "cumulative" selection. Therefore, for 
example, for each cluster the 10 most luminous red-sequence 
galaxies are present in all the other samples with larger num- 
ber of galaxies. On the other hand, when a cluster field is 
spectroscopically observed, completeness to a given limiting 
magnitude is not always achieved. In fact, the number of slits 
per mask in multi-slit spectrographs is fixed, hence the cen- 
tral, high-density regions of galaxy clusters can often be sam- 



pled only to brighter magnitudes than the external regions. 
As a consequence, the spatial distribution of the galaxies se- 
lected for spectroscopy turns out to be more extended than the 
parent spatial distribution of cluster galaxies. In the analy- 
ses presented here, we do not model this observational li mita- 
tion. I ndeed, as described in the companion paper Bazin et al.] 
(I2012h . such difficulty could be easily overcome by apply- 
ing multiple masks to the same field, which would allow one 
to achieve high completeness. For each cluster and for all 
the three orthogonal projections, we then compute the robust 
estimation of the velocity dispersion Beers et al. (1990) with 
different numbers of galaxies and compare it with the intrinsic 
ID velocity dispersion. Fig. [8] shows the probability distribu- 
tion function (PDF) of the ratio between the velocity disper- 
sion computed with different numbers of bright red- sequence 
cluster galaxies (dNgai) and the intrinsic ID velocity disper- 
sion ((Jid) obtained by stacking results from all the lines of 
sight of the cluster sample. Different colors refer to different 
numbers of galaxies and the mean of each distribution is high- 
lighted at the top of the plot with a vertical line segment. We 
note that when large numbers of galaxies are used to estimate 
the velocity dispersion, the probability distribution function 
is well represented by a log-normal distribution centered at 
zero. As a result dynamical masses obtained from large num- 
bers of bright red-sequence cluster galaxies are unbiased with 
respect to the intrinsic ID dynamical mass. However, when 
the number of red- sequence galaxies used to estimate the ve- 
locity dispersion is lower than ~ 30, the corresponding PDF 
starts to deviate from a symmetric distribution and its mean 
is biased towards lower values. This effect is evidence of a 
dynamically cold population of luminous red galaxies whose 
velocities are significantly slowed due to dynamical friction 
processes. Indeed dynamical friction is more efficient for 
more massive galaxies, hence the velocity bias is expected 
to be more important for the bright end of the ga l axy popula- 
tion ( e.R. Biviano et al. 1992, Adami e t 10119981 iCappiet al.l 
120031 lGoto..2005, .Biviano et al. 2006). 

To verify this we compute aNgai in the same way described 
above, but starting from galaxies that are randomly selected 
with respect to luminosity. Note that in this case we only 
randomly select galaxies, but we don't change the "cumula- 
tive nature" of our selection and the subsequent estimated ve- 
locity dispersion when using larger numbers of galaxies. We 
then calculate the corresponding dynamical masses in the case 
of galaxies selected according to the procedure described in 
Sect. Eland in the case of random selection using the different 
number of galaxies listed in Table [3] The resulting stacked 
dynamical masses for the full sample of clusters and for the 
three orthogonal projections are shown in Fig. |9]as a function 
of the intrinsic virial mass Mvir • The dashed purple-black 
line is the one-to-one relation and the solid green lines are the 
median and 16 and 84 percentiles of the distributions. The 
left panel of Fig. [9] represents the original distribution, while 
the right panel represents the randomly selected distribution. 
As expected from Fig.[8l if velocity dispersions are computed 
from red- sequence galaxies selected according to their lumi- 
nosity, a clear bias is introduced in the estimated dynamical 
mass. Moreover we can see that the distance between the me- 
dian line and the 84 percentile line is smaller than the distance 
between the median and the 16 percentile line, because the 
distribution is no longer symmetric. 

Furthermore it appears that the bias present in the estimated 
dynamical mass does not depend on the cluster mass. On the 
other hand, if we randomly select galaxies (right panel), the 



Saro et al. 













iiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiii 



iiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiii 



iiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiii 




-22 -21 -20 -19 -18 -17 -22 -21 -20yQq-19 -18 -17 -22 -21 -20 -19 -18 -17 

Fig. 7. — Color magnitude relation for all the galaxies within i^vir^t six different redshifts. Color-magnitude relations are expressed as ^f — r as a function of 
g for redshift z< 0.35, r — i as a function of i for redshifts 0.35 <z< 0.75 and as i — z as a function of z for redshifts larger than 0.75 (see text for further 
details). Symbols with different colors refers to different galaxy clusters in each separate redshift bin. Dashed black-red line is the model given by Song et al.i 
|2011). The solid green lines are the best fit to the simulated red-sequence relation used in this work and dashed green lines enclose 68% of the galaxies. The 
area between them represents the color space used for the selection of galaxies described in Sect. |4] 
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Fig. 8. — The probability distribution function of the measured velocity 
dispersion computed with different numbers of red-sequence cluster galaxies 
sorted by luminosity and normalized by the intrinsic ID velocity dispersion 
using the full galaxy sample. The position of the mean of each curve is high- 
lighted with a vertical line segment at the top of the figure. Small samples of 
the most luminous galaxies exhibit biases in dispersion and significant asym- 
metries in the PDF. 



bias is reduced, and we obtain a more symmetric distribu- 
tion. We also check for a possible redshift dependency on 
the velocity bias or dynamical friction. For this purpose we 
split our sample of clusters into two different redshift bins 
and show separately in Fig. [TOlthe relation between the true 
cluster virial mass and the estimated dynamical masses com- 
puted with different number of bright red-sequence galaxies 
selected according to their luminosity in the case of low red- 
shift (left panel) and high redshift (right panel) clusters. Ob- 
viously the number of clusters and their mass distribution is 
a strong function of redshift. However, it is worth noting that 
the impact of dynamical friction on the estimation of veloc- 
ity dispersion and dynamical mass does not vary much with 
cluster mass or redshift. 

Using the results of these mock observations we express 
both the velocity dispersion bias, represented by the position 
of the vertical segment at the top of Fig. [H and the charac- 
teristic width of each distribution shown in Fig. [8] with the 
following parametrisation: 

< ln{aNgai/(TiD) >= 0.05 - O.bl/y/Ngal, (6) 



^ln{(TNgal/criD) 



-0.037+ 1.047/ yiV^. (7) 
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Fig. 9. — The relation between Mvir and the dynamical mass for all the clusters in the sample and for each orthogonal projection. For each cluster the 
dynamical mass is infer-ed by applying Eq. |2]to the robust estimation of the velocity dispersion computed using different number of galaxies (Table [3). Left 
panel is for bright red-sequence galaxies sorted according to their luminosity and right panel is for a randomly sorted array of galaxies. Dashed purple-black line 
is the one-to-one relation, and solid green lines represent the 16, 50 and 84 percentiles. 



This parametrisation is valid only in the limit of the number of 
galaxies used in this study (between 10 and 100). For example 
if the dynamical mass is estimated starting from a number of 
galaxies larger than 100, the bias would presumably be zero 
rather than negative as implied from Eq. [6] 

We demonstrate in Fig. [TT] that by applying Eq. [6] to the 
velocity dispersion estimated with different numbers of red- 
sequence cluster galaxies we are able to remove the bias 
induced by the dynamical friction. In particular, Fig. [TT] 
shows the relation between true virial mass and the dynam- 
ical mass estimated using the most luminous 100, 50 and 15 
red-sequence galaxies. Dynamical masses are computed by 
applying Eq.[2l directly to the velocity dispersion (left panels) 
and to the velocity dispersion corrected according to Eq. |6] 
(right panels). Dynamical friction is affecting mostly the 
bright end of the red-sequence cluster galaxies population and 
therefore the bias is larger in the case of the smallest sample 
(lower left panel). Consistently the correction given by Eq.[6] 
is larger in this case, whereas it is negligible in the other cases 
(50 and 100). 

4.2. Impact ofPoisson Noise 

In this work we restrict our analyses to all the galaxies 
with stellar masses predicted by the adopted SAM larger than 
5 X 10^ M0. The total number of galaxies within the virial 
radius i^vir is therefore quite large and even for the poorer 
clusters with Mvir ^ 10^^ M©, the number of galaxies used 
to compute the ID velocity dispersion is of the order of 
NiD ^ 200. As a result, in the absence of any dynamical 
friction effect, the associated characteristic scatter to the ra- 
tio cFNgai/cTiD is wcll represented by the Poissonian factor 
y^2Ngai. To demonstrate it, we show in Fig.[T2]the evolution 
of the scatter in the relation between the true virial masses 
and the dynamical masses as a function of redshift. For each 
cluster, dynamical mass is estimated starting from the veloc- 



ity dispersion of the 100 most luminous red-sequence galax- 
ies through Eq. O The resulting scatter is highlighted as a 
cyan solid line. We also show the evolution of associated scat- 
ter when dynamical mass is computed from the intrinsic 3D 
(dashed red line) and ID (dotted black line) velocity disper- 
sions. Moreover, similarly to Fig. |4] we separately show the 
predicted scatter obtained by adding in quadrature the scatter 
associated to the ID velocity dispersion with the Poisson term 
^J2Ngal (dashed-triple dotted green line) or with the factor 
given by Eq. [7] (dashed-dotted purple line). We note that, as 
expected, both predictions agree very well with the measured 
evolution of the scatter. 

However, if a lower number of galaxies is used to calcu- 
late the dynamical mass, a difference in the two predictions 
emerge. For example, in Fig. [13] we show the same compu- 
tation highlighted in Fig. [121 but with a number of galaxies 
equal to 50 (left panel) and 15 (right panel). We note in par- 
ticular that the observed evolution of the scatter of the rela- 
tion among the virial mass and the dynamical mass is well 
described by adding in quadrature to scatter associated to the 
intrinsic ID dynamical mass the term given by the fitting for- 
mula of Eq. [7] On the contrary, if only the Poisson term 
y^2Ngai is taken into account, the predicted scatter is un- 
derestimated with respect to the measured one. Furthermore, 
note that while on the one hand in Figures [TT]and[8] we showed 
that the dispersion calculated using 50 galaxies does not in- 
troduce a significant mass bias; however, it is clear that the 
Poisson term is no longer adequate to explain the real scat- 
ter. On the other hand, if the contribution to the scatter due 
to dynamical friction is included through Eq. [7] we are able to 
recover the right answer. 

4.3. Impact of Interlopers 

Finally, to have a more coherent and realistic approach to 
our analyses, we further investigate the effect of interlopers as 
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Fig. 10. — Same as for the left panel of Fig. |9] but dividing our cluster sample in 2 redshift bins. Left panel is for z < 0.5 and right panel is for z > 0.5. 



a possible source of systematics in the computation of clus- 
ters dynamical mass. For this purpose, for each snapshot 
and projection, we construct a number of cylindrical light- 
cones centred at each cluster with height equal to the full 
simulated box-length and different radius spanning the inter- 
val 0.2 to 2.4 i^vir- The different aperture values used are 
listed in Tabled We then apply an initial cut of 4000 km s~^ 
to select galaxies within the cylinders. For each cylindrical 
light-cone realisation, we then initially select a different num- 
ber of members in the color-magnitude space described in 
Sect. |4] ranging from 10 to 100 galaxies as shown in Table 
[51 Several techniques have been developed to identify and 
reject interlopers. Such methods have been studied before 
typically using randoml y selected d ark m atter particles (e.g . 
IPerea et al. 1990; Diafe rio & Gelledll997l: iLokas et al.ir2006l; 
[Woitak et al. 2007, 2009^ and more recently using subhalos 
by Bivianoetal. (2006) and White et al. (2010). However, 
for the purpose of this work, here we simply apply a 3 a clip- 
ping procedure to its robust estimation of the velocity disper- 
sion [B eers et al. (1990) to reject interlopers, as discussed in 
iBazin et a l. (2012). This leads to a final spectroscopic sample 
of galaxies for each cluster, at each redshift, for each pro- 
jection, within each different aperture and for each different 
initially selected number of red- sequence galaxies. Fig. [14] is 
a schematic representation of the procedure we follow to ob- 
tain from each cluster and projection different estimation of 
the velocity dispersion according to different "observational 
choices". 

From final spectroscopic sample of galaxies described 
above, we compute the fraction of interlopers (arbitrary de- 
fined here as galaxies lying at a cluster centric distance larger 
than 3 x i^vir) as a function of the aperture by stacking to- 
gether the sample in different bins according to their redshift, 
to the number of galaxies used to evaluate their velocity dis- 
persion and to the cluster masses. This can be seen in Fig- 
ure [T5l and is in go od agreement with previous works (e.g. 
iMamon et al.ll2010h . The two upper panels and the lower-left 



panel show the fraction of interlopers as a function of aperture 
respectively color coded according to the number of galaxies 
(panel A), to the redshift (panel B) and to the cluster mass 
(panel C). As expected, the fraction of interlopers rises with 
the aperture within which the simulated red- sequence galaxies 
were initially chosen. This indicates that even red sequence, 
spectroscopically selected samples are significantly contami- 
nated by galaxies lying at distances more than three times the 
virial radius from the cluster. 

On the other hand a much weaker dependency between the 
number of selected red sequence galaxies and the fraction of 
interlopers is highlighted on the upper-left panel of Fig. [15] 
(A). Whether one has small or large samples the fraction of 
interlopers remains almost the same. The upper-right and the 
bottom-left panels are showing that the fraction of interlop- 
ers is larger at larger redshifts consistently with a denser Uni- 
verse (B), and is a steeper function of aperture for lower mass 
clusters (C). Since in the hierarchical scenario more massive 
halos forms at later times than the lower mass ones, these two 
variables are clearly correlated. Thus, we also show in the 
bottom-right panel labeled "D" how the fraction of interlop- 
ers varies as a function of redshift by stacking together the 
sample in different mass bins. Most massive clusters are not 
formed yet at high redshifts, therefore above certain redshifts 
the redder lines go to zero. Although oscillating, an evident 
tendency of increasing fraction of interlopers is associated to 
larger redshift, whereas at fixed z there is no clear dependency 
of the fraction of interlopers from the clusters mass. We stress 
however that all the relations shown in Fig. [15] are meant to 
describe the qualitative dependency of the interloper fraction 
from the analysed quantities. 

In a similar way we compute the mean velocity bias 
(defined as the ratio between the measured velocity dis- 
persion and the intrinsic line-of-sight velocity dispersion: 
^{N auR±Mvir,z)/^iD) as a function of the aperture by stack- 
ing together the sample in different bins according to their 
redshift, to the number of spectroscopic galaxies and to the 
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Fig. 11. — The relation between Mvir and the dynamical mass for all the clusters in the sample and for each orthogonal projection. Left panels: For each 
cluster the dynamical mass is infer-ed by applying Eq. |2]to the robust estimation of the velocity dispersion computed using the 100, 50 and 15 most luminous 
red-sequence cluster galaxies (with distance from the centre smaller than i^vir) and show respectively in the upper, middle and lower panels. Right panels: Same 
as for left panels, but velocity dispersions are corrected according to Eq. |6] Dashed purple-black line is the one-to-one relation, and solid green lines represent 
the 84, 50 and 16 percentiles. 
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Fig. 12. — The evolution of the 1 a scatter as a function of redshift in log- 
space for the following quantities. Dashed red (dotted black) line is for the 
ratio between the estimated dynamical mass Mn^j^, {Mir)) computed from 
the 3D (ID) velocity dispersion and the virial mass Mvir • The solid cyan 
line is for the ratio between the measured dynamical mass M^yn computed 
from the 100 most luminous red-sequence galaxies within the virial radius of 
each cluster along each line of sight and the virial mass Mvir • The dashed- 
dotted green line is the expected scatter in mass obtained by multiplying B 
by the term given by adding in quadrature the scatter from the ID velocity 
dispersion and a Poissonian term equal to V2 x 100. The dashed-dotted 
purple line is the expected scatter in mass obtained by multiplying B by the 
term given by adding in quadrature the scatter from the ID velocity dispersion 
and a term computed with the fitting formula of Eq.|7] 



cluster mass. This can be seen in Figure [161 The two up- 
per panels and the lower-left panel show the velocity bias as a 
function of aperture respectively color coded according to the 
number of galaxies (panel A), to the redshift (panel B) and to 
the cluster mass (panel C). Interestingly, the velocity bias has 
a minimum when velocity dispersions are evaluated within 
R± ~ i^vir, and rises at both smaller and larger radii. In par- 
ticular, for projected radii < i^vir, where the effect of interlop- 
ers is smaller, we recover the expe cted decrease of th e average 
velocity dispersion profile (e.g. Bivia no et al.ll2006h as a func- 
tion of aperture. On the other hand, for R± > i^vir, the larger 
contamination from interlopers is significantly affecting and 
boosting the velocity bias. Furthermore, as expected, dynam- 
ical friction is also affecting the estimated velocity dispersion, 
when the latest is computed with a small number of selected 
red sequence galaxies (A). Indeed, by applying Eq.[6]to the es- 
timated velocity dispersion we are able to successfully remove 
the degeneracy between the velocity bias and the number of 
galaxies within projected aperture R± < i^vir (Fig- El- The 
upper-right and the bottom-left panels of Fig. [16] show that, 
consistent with the fraction of interlopers, the velocity bias 
computed within R± > i^vir is larger at larger redshifts (B) 
and is a steeper function of aperture for lower mass clusters 
(C). Finally, a mild dependence of redshift for fixed mass is 
highlighted in the bottom-right panel (D). 

To better understand how interlopers affect the inferred ve- 
locity dispersion we select as an example all clusters with 
Mvir larger than 5 x 10^^ M©. For each of the three or- 
thogonal projections we then initially select the most lumi- 
nous 25 red-sequence galaxies as described in Sect. |4] within 
a projected distance of l.Si^vir- We then apply the same pro- 
cedure described above to reject interlopers and obtain a final 
list of galaxies. From this list of galaxies we then identify 



the "true" cluster members and the interlopers. We show in 
the left panel of Fig. [18] a map representing the stacked dis- 
tribution of the velocity of the cluster galaxies as a function 
of the projected separation from the cluster center R±/Ryir. 
Note the typical trumpet shape of the e xpected caust ic distri- 
bution (Diaferio & Geller 1997l lSerrae t al. 201 1 , Zha ng et al.] 
1201 lb . On the top of this map, we overplot as contours the 
stacked distribution of the interloper population that the So- 
clipping procedure was not able to properly reject. A large 
fraction of high velocity interlopers are still present after fore- 
ground and background removal and thus they will bias high 
the estimated velocity dispersion. 

This map highlights how caustic based techniques are po- 
tentially more effective to remove interlopers than a simple 
3(7 clipping. However, observationally, a much larger number 
of galaxies than the 25 spectra used here is typically needed 
to apply these more sophisticated methods. 

We also show in the right panel of Fig. [18] as solid black 
and dashed red histograms respectively the distribution of ve- 
locities for both the cluster galaxies and the interlopers pop- 
ulation. The expected Gaussian velocity distribution is over- 
plotted as a solid black Gaussian with a standard deviation 
given by Eq. [6] and Ngai = 25. The absolute normalisa- 
tions of the histograms are arbitrary, but the relative ratio of 
the two histograms is representative of the ratio between the 
number of cluster galaxies and interlopers. Note also that a 
large fraction of low velocity interlopers is present. These in- 
terlopers are mostly red-sequence galaxies which lie at about 
the turn-around radius of the cluster over-density and there- 
fore have associated redshifts which are consistent with the 
cluster redshift. As discussed above, a simple 3a clipping 
technique is not able to effectively remove high velocity in- 
terlopers, and therefore is biasing high the inferred velocity 
dispersion. On the contrary caustic based methods are able to 
remove this high velocity interlopers population, but are not 
effective to reject this low velocity galaxies at around the turn- 
around radius. As a net result, velocity dispersions computed 
after interlopers rejection based upon caustic techniq ues will 
be biased low ( Woitak et al.ll2007LIZhang etaniMTTl) . 

As mentioned above, for each cluster along all the pro- 
jections we end up with different samples of red- sequence 
galaxies that the 3a clipping procedure recognises as "spec- 
troscopic members". Therefore, for each different initially 
selected number of red- sequence galaxies, we measure the 
robust estimation of the velocity dispersion. We then ap- 
ply Eq. [2] to estimate the dynamical mass. Left panel of 
Fig. [19] shows the corresponding relation between the result- 
ing dynamical mass and the true virial mass for all the sample 
stacked together. The dashed black-purple line is the one to 
one relation, whereas the green lines show the 16, 50 and 84 
percentiles. Note that the sample shown here is volume lim- 
ited, and so the distribution in mass is different than the typi- 
cal observational samples. Furthermore, the same clusters ap- 
pear several times with dynamical masses computed from dif- 
ferent number of galaxies on each projection, and within dif- 
ferent projected radii at all the redshifts. When red-sequence 
galaxies are selected within a projected radius from a light- 
cone regardless of their true 3D distance from the centre of the 
cluster, the relation between the virial mass and the inferred 
dynamical mass is much broader. In particular, by looking at 
the median of the distribution, it is possible to notice that a 
systematical overestimation of the dynamical mass is present 
at all cluster masses, as expected from the interloper contribu- 
tion previously discussed. Furthermore, especially at the low 
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Fig. 13. — Same as for Fig. [12] but with a number of galaxies respectively equal to 50 (left panel) and 15 (right panel). 









a 


cr' 


y^ ~~^ 


\l 


Ngal=25 


o 


/ a 


\ 


r<0.8R. 


.£2 


Ngal=25 


) 




c 


Vr<0.6R 


/ 


• . 


CJ 


\ " 


% 




.c 


^^^S^^^^^y^ 




A 


1/ 


o \ 




f NgaA 
I 25 i 




Ngal=25 


" 




Q. \ 


<0.4R. / 








/ 


a 


X 




Ngal=25 ) 



\ 



% 



Cluster 

\ z, M .projection/ 



^ Ngal V 



/ Ngal \ 

)X/ 

X a \ 





Ngal=15 , 

r < 0.2 R^. / 


Ngal=15 
\ r<0.4R. / 


^^^^ ^ 



/ a 

Ngal=15 
\ r<0.6R. 



Fig. 14. — Schematic representation of the parameter space explored in 
this work. See Table [3]f or specific ranges in each parameter. 



mass end of the cluster galaxy distribution, the presence of a 
significant population of catastrophic outliers is making the 
relation among virial mass and dynamical mass very asym- 
metric and causing a severe boosting of the dynamical mass. 
These outliers are likely related to cases where the sim- 
ple 3cr clipping procedure is not sophisticated enough to ef- 
fectively separate the foreground and background interloper 
galaxies from the proper cluster galaxies. To verify this hy- 
pothesis we show in the right panel of Fig. [19] the same com- 
putation as for the left panel, but restricting our sample to 
only the cases in which the presence of interlopers is smaller 
than 5%. We note how this sub-sample qualitatively looks 
very similar to left panel of Fig. [9] which by construction 
contains only cluster galaxies. Furthermore, once the con- 
tribution from interlopers is removed, the bias of dynamical 
mass over the true mass disappears. However, remember that 
Fig. [19] shows that without interlopers dynamical masses are 
on average underestimated compared to the true virial mass, 
as expected from the effect of dynamical friction described 
in Sect. 14.11 Moreover, as expected from the lower panels of 



Fig. \T5\ the adopted interlopers rejection method is more ef- 
fective for more massive clusters. Clearly the interloper effect 
on the dynamical mass is more severe at the low mass end of 
the cluster population. 

Because the color selection of cluster members is a crucial 
point in this analysis, the results presented here obviously de- 
pend on the adopted galaxy formation model at some level. 
On the one hand it is true that the model is not perfectly re- 
producing the observed properties of the cluster galaxy popu- 
lation. On the other hand we also do not take into account any 
observational uncertainty which will instead affect the real 
data, for example broadening the observed red-sequence at 
fainter magnitudes. 

To estimate the sensitivity of the color selection to uncer- 
tainties in the galaxy modeling on the above described results, 
we select red- sequence galaxies with a different criteria than 
the one described in Sect. |4] Instead of selecting the area in 
color-magnitude space which encompasses 68% of the clus- 
ter galaxies, we select all galaxies within a fixed ±0.15 mag 
along the fitted red- sequence re lation, simila rly to the adopted 
criteria in the companion paperlBazin et al.lOOll ). This is on 
average a factor of ~ 2 in magnitude larger than the former 
threshold (depending on the redshift ranging from ~ 0.5 — 3, 
as highlighted in Tab. [2]). Then, we reject interlopers and com- 
pute velocity dispersions and subsequent dynamical masses as 
described in the above sections. We find that the fraction of 
interlopers which the 3a clipping procedure is not able to re- 
ject is on average in agreement within ~ 3% with the previous 
color selection. In particular for clusters with Mvir larger than 
4 X lO^^Mvir the agreement is better than 1%. 

We show in the left panel of Fig.[T9]the resulting 16, 50 and 
84 percentiles overplotted as red continuous lines. We note 
that a larger effect from the interlopers is present in compari- 
son with the previous analyses, as expected from the broader 
color selection adopted. In particular, larger differences ap- 
pear at the low mass end of the cluster galaxy population, 
where a significant increase of catastrophic outliers in the 
overestimation of the dynamical mass is visible. On the other 
hand, the average population is not affected by much. As a 
net result, a changing in the color selection of a factor ~ 2 
implies a change in the estimated velocity dispersion by less 
than ^ 3%. In particular, this difference reduces to less than 
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panel labeled A), redshift (top right panel labeled B) and mass of the cluster (bottom left panel labeled C). Bottom right panel: The stacked mean fraction of 
interlopers as a function of redshift color coded according to the cluster mass (labeled D). 




1% for clusters with Mvir larger than 5 x 10^"^ M(^ 
4.4. Unbiased Dispersion Mass Estimator 



Similarly to Sect. 14.11 we try to parametrise as a function of 
the variables in Table [3] (aperture, number of spectra, redshift 
and cluster mass), the way that interlopers affect the inferred 
dynamical mass. However, we could not find a satisfactory 
analytical solution to easily model the measured velocity dis- 
persion of clusters as a function of the above described vari- 
ables, due to the non-linear interplay of the explored param- 
eter space highlighted in Fig. [161 Therefore, we numerically 
compute the mean and the associated standard deviation of 
the ratio between the observed and the ID intrinsic velocity 
dispersion in different bins of the parameter space as high- 
lighted in TableO In this way, given the cluster mass, redshift 
and the number of red- sequence galaxy spectra within a given 
projected radius used to compute the velocity dispersion, we 
can correct for the average bias affecting the estimation of the 
dynamical mass. We show in Fig. [20lthe same relation de- 
scribed in the left panel of Fig. [19] when such corrections are 



included. We remark that the bias is effectively removed at 
all the mass scales analysed here. Furthermore, by compar- 
ing the 84 percentile and median lines at the low mass end of 
the left panel of Fig. [19] with the ones in Fig.[20l we note that 
while the former are separated by about an order of magnitude 
in dynamical mass, for the later this difference is reduced to 
about 0.8 dex. 

5. DISCUSSION AND CONCLUSIONS 

We have examined the use of velocity dispersions for un- 
biased mass estimation in galaxy clusters using the publicly 
available galaxy catalogue produced with the semi-analytic 
model by De Lucia & Blaizot (2007) coupled with the N-body 
cosmological Millennium Simulation (Springel et al. 2005). 
In particular, we selected all galaxies in the SAM with stellar 
mass larger than 10^ Mq and analysed a sample consisting of 
more than -- 20000 galaxy clusters with Mvir > 10^"^ M© up 
toz- 1.2 (Tab: [B. 

First we explore the properties of the full galaxy sample and 
then we increase the level of complications to mimic the spec- 
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troscopic selection that is typically undertaken in real world 
studies of clusters. Then we work through a series of con- 
trolled studies in an attempt to disentangle the different ef- 
fects leading to biases and enhanced scatter in velocity dis- 
persion mass estimates. Ultimately our goal is to inform the 
dispersion based mass calibration of the SPT cluster sample 
iBazin et al. 2012), but we explore a broad range in selection 
in hopes that our results will be of general use to the commu- 
nity. 

Our primary conclusions for the full subhalo population 
are: 



• We measure the galaxy (i.e. subhalo) velocity dis- 
persion mass relation and show that it has low scatter 
(~ 0.14 in /n(M)) and that subhalo dispersi ons are 
<3% lower than DM dispersion in Evra rd et al.l (|2008). 
This difference corresponds to a < 10% bias in mass for 
our halos if the DM dispersion- mass relation is used, 
and is consistent with previous determination of sub- 
halo velocity bias. 



• We explore line of sight velocity dispersions of the 
full galaxy populations within the cluster ensemble and 
confirm that the triaxiality of the velocity dispersion el- 
lipsoid is the dominant contributor to the characteris- 
tic ~35% scatter in dispersion based mass estimates. 
We show that this scatter increases with redshift as 
(7{z) c^0.3 + 0.075z. 

• We measure the principal axes and axial ratios of the 
spatial galaxy distribution ellipsoid, showing that there 
is a slight (~ 5%) preference for prolate distributions; 
this property has no clear variation with mass or red- 
shift. We examine the line of sight velocity dispersions 
along the principle axes, showing that the slight prefer- 
ence toward prolate geometries translates into a slight 
(~ 1%) bias in the dispersion mass estimates extracted 
from line of sight measures. 

Our primary conclusions for the spectroscopic subsamples 
of subhalos are: 
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Fig. 17. — Upper panels and bottom left panel: Same as Fig. [16]- panel 
A, but corrected for dynamical friction according to Eq. |6] 

• We characterize the bias (Eqn. [6]) and the scatter 
(Eqn. [7]) in the Hne of sight velocity dispersion intro- 
duced by selecting a subset Ngai of the most luminous 
red sequence galaxies within a cluster. The bias is sig- 
nificant for samples with Ngai < 30 and is likely due 
to dynamical friction of these most massive subhalos. 
The scatter cannot be fully explained through a combi- 
nation of intrinsic scatter in the relation between mass 
and the 3D dispersion of all galaxies (i.e. departures 
from equilibrium), scatter of the line of sight dispersion 
around the 3D dispersion (halo triaxility) and Poisson 
noise associated with the number of subhalos Ngai . A 
further component of scatter due to the presence of a 
dynamically cold population of luminous red- sequence 
galaxies is needed to explain the full measured scatter. 

• We explore the impact of interlopers by creating spec- 
troscopic samples using (1) red sequence color selec- 
tion, (2) a maximum projected separation from the clus- 
ter center, and (3) A^-sigma outlier rejection in line of 
sight velocity. In these samples the interloper frac- 
tion (contamination) can be significant, growing from 
~ 10% at the projected virial radius to ~35% at twice 
the project virial radius. The contamination fraction has 
a much weaker dependency on the sample size Ngai. 
We explore the dependence on mass and cluster red- 
shift, showing that within a fixed aperture, contamina- 
tion is a factor of '^ 2 worse at redshift z ^ 1 than at 
z = 0. Furthermore, we show that the fraction of in- 
terlopers is a steeper function of aperture for low mass 
clusters, but that at fixed redshift contamination does 
not change significantly with mass. We show that con- 
tamination is significant even if a more sophisticated 
caustic approach is used to reject interlopers, demon- 
strating that even clusters with large numbers of spec- 
troscopic redshifts for red sequence selected galaxies 
suffer contamination from non-cluster galaxies. We fur- 
ther study how interlopers are affecting the estimated 
velocity bias. We find that the velocity bias has a 
minimum if computed within R± ~ Rvir- This is 
due to the balancing effect of larger intrinsic velocity 
bias at smaller radii and larger contamination at larger 



radii. Furthermore, we show that if velocity dispersions 
are computed within projected aperture R± larger than 
^ i^vir, the velocity bias is a steeper function of R± 
for higher redshifts and lower cluster mass, as expected 
from the contamination fraction. 

• We study how changing the color selection affects the 
fraction of interlopers and the subsequent effect on the 
estimated velocity dispersion and dynamical masses. 
We find that doubling the width of the color selection 
window centered on the red sequence has only a mod- 
est impact on the interloper fraction. The primary ef- 
fect of changing the color selection is on the filtering of 
catastrophic outliers. This results in changes to the es- 
timated velocity dispersion virial mass relation at the 
level of 1% in mass. We also show that uncertain- 
ties in the color selection are more important for low 
mass clusters than for the high mass end of the clus- 
ter population, which is because the dispersions of low 
mass clusters are more sensitive to catastrophic outliers. 
The rather weak dependence of the dispersion based 
mass estimates on the details of the color selection sug- 
gests also that uncertainties in the star formation histo- 
ries (and therefore colors) of galaxy populations in and 
around clusters are not an insurmountable challenge for 
developing unbiased cluster mass estimates from veloc- 
ity dispersions. 

• We present a model to produce unbiased line of sight 
dispersion based mass estimates, correcting for inter- 
lopers and velocity bias. We also present the probabil- 
ity distribution function for the scatter of the mass esti- 
mates around the virial mass. These two data products 
can be used together with a selection model describing 
real world cluster dispersion measurements to enable 
accurate cluster mass calibration. 



In a companion paper, iBazin et aH (120 12h apply this model 
in the dispersion mass calibration of the SPT Sunyaev- 
Zel'odovich effect selected cluster sample. We identify the 
following key remaining challenges in using dispersions for 
precise and accurate mass calibration of cluster cosmology 
samples. Surprisingly, the larger systematic uncertainty has 
to be ascribed to our relative poor knowledge of the velocity 
bias between galaxies or subhalos and DM. A conservative 
estimate of this systematic is at the level of < 5% and arises 
from the comparison of different simulati ons and different al - 
gorithms for subhalo identification (e.g. lEvrard et"ani2008h . 
The systematic uncertainty in the color selection of galaxies 
and its subsequent mapping between line of sight velocity dis- 
persion and mass is at relatively smaller level. Indeed we can 
estimate it a t a < 1% le v el fo r samples selected as the ones 
described in iBazin et aH ('2012), despite the fact that galaxy 
formation models involve a range of complex physical pro- 
cesses. In other words, systematics in predicting galaxy prop- 
erties (e.g. luminosity, colors, etc.) due to subgrid physics 
associated with magnetic fields, AGN and supernova feed- 
back, radiative cooling and the details of star formation, do 
not appear to significantly change the spectroscopic sample 
selection. On the other hand, simulations including different 
physical treatments of gravity are affecting the dynamics of 
the spectroscopic selected sample at a higher level than we 
expected. Given that the current dominant contributor to the 
systematics floor is an issue associated with the treatment of 
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Fig. 18. — Left panel: The color map represents the distribution of the line of sight velocity of cluster galaxies (within 3i?vir) normalized to the intrinsic 3D 
velocity dispersion of clusters as a function of the projected distance from the cluster center in units of i?vir for the sample described in the text. The contour 
lines represent the same distribution for the interloper galaxies. Right panel: The distribution of velocities in units of the intrinsic 3D velocity dispersion for 
the cluster galaxy population (solid black histogram) and for the interloper population (dashed red histogram). The normalisation is arbitrary, while the relative 
ratio of the two histograms reflects the sample described in the text. The solid black Gaussian is the expected distribution with width given by the Eq.|6]and 
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gravitational physics, there are reasons to be optimist that fu- 
ture simulations will be able to reduce the current systematics 
floor. 



We acknowledge Jonathan Ruel for very useful discus- 
sions and support from the Deutsche Forschungsgemeinschaft 
funded Excellence Cluster Universe and the trans-regio pro- 
gram TR3 3: Dark Universe. 



REFERENCES 



Abell, G.O. 1958,ApJS,3,211 

Adami, C, Biviano, A., & Mazure, A. 1998, A&A, 331, 439 
Andersson, K., Benson, B. A., Ade, R A. R., et al. 2011, ApJ, 738, 48 
Barnes, J., & Hut, R 1986, Nature, 324, 446 



Bazin et al. 2012, In Preparation 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Benson, A. J. 2010, Phys. Rep., 495, 33 



18 



Saro et al. 



Full sample 



with corrections 

T 




14.6 14J 

log Kw [ M© ] 

Fig. 20. — Same as for the left panel of Fig. [19] but velocity dispersions 
were numerically corrected as described in the text. Dashed purple-black 
line is showing the oat-to-one relation, while solid green lines are the 16, 50 
and 84 percentile. Note that with these corrections the dynamical mass is an 
unbiased estimator of the true mass. 

Biviano, A., Girardi, M., Giuricin, G., Mardirossian, R, & Mezzetti, M. 

1992, ApJ, 396, 35 
Biviano, A., Murante, G., Borgani, S., et al. 2006, A&A, 456, 23 
Bohringer, H., Voges, W., Huchra, J. P., et al. 2000, ApJS, 129, 435 
Cappi, A., Benoist, C, da Costa, L. N., & Maurogordato, S. 2003, A&A, 

408, 905 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 

371 
De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 

De Lucia, G., Kauffmann, G., Springel, V., et al. 2004a, MNRAS, 348, 333 
De Lucia, G., Kauffmann, G., & White, S. D. M. 2004b, MNRAS, 349, 1101 
Diaferio, A., & Geller, M. J. 1997, ApJ, 481, 633 
Dressier, A., & Shectman, S. A. 1988, AJ, 95, 985 
Evrard, A. E., Bialek, J., Busha, M., et al. 2008, ApJ, 672, 122 
Faltenbacher, A., & Diemand, J. 2006, MNRAS, 369, 1698 
Fassbender, R., Bohringer, H., Nastasi, A., et al. 2011, New Journal of 

Physics, 13, 125014 
Gao, L., De Lucia, G., White, S. D. M., & Jenkins, A. 2004, MNRAS, 352, 

LI 



Geller, M. J., & Beers, T. C. 1982, PASP, 94, 421 

Gladders, M. D., & Yee, H. K. C. 2005, ApJS, 157, 1 

Goto, T. 2005, MNRAS, 359, 1415 

Hao, J., McKay, T. A., Koester, B. R, et al. 2010, ApJS, 191, 254 

Jing, Y. R, & Suto, Y. 2002, ApJ, 574, 538 

Kasun, S. R, & Evrard, A. E. 2005, ApJ, 629, 781 

Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 1999, 

MNRAS, 303, 188 
Koester, B. R, McKay, T. A., Annis, J., et al. 2007, ApJ, 660, 221 
Lokas, E. L., Wojtak, R., Gottlober, S., Mamon, G. A., & Prada, R 2006, 

MNRAS, 367, 1463 
Mamon, G. A., Biviano, A., & Murante, G. 2010, ArXiv e-prints 
Marriage, T. A., Acquaviva, V., Ade, R A. R., et al. 2011, ApJ, 737, 61 
Mohr, J. J., Evrard, A. E., Fabricant, D. G., & Geller, M. J. 1995, ApJ, 447, 

8-F 
Pacaud, R, Pierre, M., Adami, C., et al. 2007, MNRAS, 382, 1289 
Perea, J., del Olmo, A., & Moles, M. 1990, A&A, 237, 319 
Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2011, ArXiv 

e-prints 
Reichardt et al. 2012, In Preparation 
Reiprich, T. H., & Bohringer, H. 2002, ApJ, 567, 716 
Serra, A. L., Diaferio, A., Murante, G., & Borgani, S. 2011, MNRAS, 412, 

800 
Song, J., Mohr, J. J., Barkhouse, W. A., Warren, M. S., & Rude, C. 2011, 

ArXiv e-prints 
Springel, V. 2005, MNRAS, 364, 1105 
Springel, V., White, S. D. M., Tormen, G., & Kauffmann, G. 2001a, 

MNRAS, 328, 726 
Springel, V., Yoshida, N., & White, S. D. M. 2001b, Nature, 6, 79 
Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 
Stanford, S. A., Eisenhardt, R R., Brodwin, M., et al. 2005, ApJ, 634, L129 

Staniszewski, Z., Ade, R A. R., Aird, K. A., et al. 2009, ApJ, 701, 32 
Sunyaev, R. A., & Zel'dovich, Y. B. 1972, Comments on Astrophysics and 

Space Physics, 4, 173 
Szabo, T., PierpaoH, E., Dong, R, Pipino, A., & Gunn, J. 2011, ApJ, 736, 21 
Thomas, R A., & Couchman, H. M. R 1992, MNRAS, 257, 11 
Tormen, G. 1997, MNRAS, 290, 411 

Suhada, R., Song, J., Bohringer, H., et al. 2012, A&A, 537, A39 
Vanderlinde, K., Crawford, T. M., de Haan, T., et al. 2010, ApJ, 722, 1180 
Vikhlinin, A., Burenin, R. A., Ebeling, H., et al. 2009, ApJ, 692, 1033 
Warren, M. S., Quinn, R J., Salmon, J. K., & Zurek, W. H. 1992, ApJ, 399, 

405 
White, M., Cohn, J. D., & Smit, R. 2010, MNRAS, 408, 1818 
Williamson, R., Benson, B. A., High, R W., et al. 2011, ApJ, 738, 139 
Wojtak, R., Lokas, E. L., Mamon, G. A., & Gottlober, S. 2009, MNRAS, 

399, 812 
Wojtak, R., Lokas, E. L., Mamon, G. A., et al. 2007, A&A, 466, 437 
Zhang, Y.-Y, Andemach, H., Caretta, C. A., et al. 2011, A&A, 526, A105 
Zwicky, R, Herzog, E., & Wild, P. 1968, Catalogue of galaxies and of 

clusters of galaxies (California Institute of Technology, Pasadena) 



